SMI-20/96 



Improved Algorithm for Bosonized Fermion 

Determinant. 

^ ; A.A.Slavnov * 

Steklov Mathematical Institute, Russian Academy of Sciences, 
Vavilov st.42, GSP-1, 117966, Moscow, Russia 



OS 



X 



February 1, 2008 



Abstract 



> 

| Following the line of ref. |TJ we propose an improved algorithm which 

allows to calculate a .D-dimensional fermion determinant integrating the ex- 
ponent of D + 1-dimensional Hermitean bosonic effective action. For a finite 

| extra dimension the corrections decrease exponentially. 

On 
+-> 

ri3 : 1 Introduction 

a: 

Recently we proposed an algorithm which allows to calculate D-dimensional fermion 
determinants by integrating a D + 1-dimensional bosonic effective action Q. The 
following representation for the square of the covariant Dirac operator was obtained: 



det (D + m 



,2 



r N 

lim exp{a% £ ^^(^(a^Or) + ft^^x) - 2<j>*J n )+ 



n=-JV+l s 



+i[^i(»)75p + m)0 n (z) - 0:(x) 75 (,D + m)0„ +1 (x)][6]- 1 + 

1 z 
+ 2^( J ° 2 ~ m2 )^n + -j={<f)* n {x)x{x) + x*(x)(f) n (x))]}d(j)* n d(j) n dx*dx- 

Here <ft n are bosonic fields defined on a 4+1 dimensional lattice which have the same 
spinorial and internal structure as the Dirac fields ip. The fifth component t to be 
denned on the one dimensional lattice of the length L with the lattice spacing b: 

L = 2Nb, -N < n < N (2) 
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The free boundary conditions in t are assumed: 

K = 0, n < -N, n > N (3) 
Lattice covariant derivative is denoted by 

D^{x) = hfj^x^ix + a„) - i){x)\ (4) 

D = 1/2 % (D; - D,) (5) 

(For simplicity we consider naive fermions, but all the construction is extended in a 
straightforward way to Wilson fermions). 

For finite L and b the eq.(|l|) should be corrected by the terms 

°^-JXD^ 2 ^) (6) 

where D are eigenvalues of the Dirac operator. 

One sees that although in the limit b — > 0, L — > oo the eq.(0) is exact, the 
convergence is not very fast. For numerical simulations it is important to make 
the corrections as small as possible, so one could use a one-dimensional lattice with 
relatively small number of sites N. Recently there were several publications, mainly 
based on the Lusher's proposal for bosonization of fermion determinants, where 
different aspects of numerical simulations in the bosonized models were discussed 
0, |H 0. Further development in ref. ||, 0, lead to algorithms with better 
convergence properties. 

In the present paper I propose a modified algorithm, using essentially the same 
idea as in ref. 0, but providing much better convergence. Instead of polynomial 
supression of finite size effects as in eq.(|6|), new algorithm provides exponential 
damping 

exp{— mL} (7) 

where m is a bare quark mass. 



2 Improved algorithm for lattice QCD 

Having in mind applications to QCD in this section we consider two fermion flavours 
interacting vectorially with the Yang- Mills fields U^. The reason to consider two de- 
generate flavours is the positivity of the square of the gauge covariant Dirac operator. 
It is convinient to present the square of fermion determinant in the followi ng form 

r 2 - 

/ exp{a 4 J2Y,Mx){D + m)^„(x)}## = det[ 75 (L> + m)] 2 = (8) 

J n=l x 

= / expja 4 ip(x)(D 2 — m 2 )ip(x)}dipdip 
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We again introduce five dimensional bosonic fields <fr(x, t) with the same spinorial 
and internal structure as ip(x). The notations above. 
We are going to prove the folowing identity 

f exp{a 4 ^^(^)P 2 -m 2 )^(x)}## = (9) 

J X 
N 

= lim exp{a% £ ^[(6- 2 )(^ +1 (a;)0„(a;) + /..c. - 2^) + 



n=-jV+l x 



— «C 



n+1 



x)-f 5 DL- 1 (2n + l)(f) n + h.c. 



+ l/2^{x)D 2 L-\2n + l) 2 b 2 <f) n 
^J^l\4rjx)x{x) + X %x)(P n ^))^bexp{-mL^b 2 n 2 }]}d4>* n d^ n dx*dx. 

7TL° 

Integration goes over the fields (f> n (x) satisfying the free boundary conditions @. 
The last term in the exponent in the r.h.s. of eq.(|9|) introduces the constraint 

^0 n nexp{-mL~Vn 2 } = (10) 

n 

Substituting the solution of this constraint to the eq.(||) one gets the representation 
for the determinant of the gauge covariant Dirac operator as the path integral of the 
exponent of the bosonic Hermitean action. Below we present a construction which 
justifies the eq.(|9[) and estimate the corrections due to the finite lattice spacing b 
and finite lattice size L. 

Our starting point is the integral 



N 
I 

n=-N+l x 

>2„2r-l 



r N 

1= exp{a% £ J2i( b ' 2 Wn + i(^ e M^Db 2 (n + l) 2 L- 1 }x (11) 

« A r i 1 <y> 



ex-p{i^ 5 Db n L }(f) n (x) + h.c. - 20*0„) + 

+i (™Yl\4,* n (x)x{x) + X \x)(P n (x))2nbexp{-mL- l b 2 n 2 }]}d(t)ld(t) n dx f dx. 
ttL° 



The operator 75D is Hermitean and it's eigenvalues are real. R.h.s. of eq. (|ll|) can 
be written in terms of eigenvectors of the operator 75 Z), which are denoted as (p a - 

N 

1= exp{b J2 Y^ ^T+^M-iD^in + lfL- 1 }^ (12) 

J n= _iv+l a 

exv{iD a b 2 n 2 L- l }<P« + h.c. - 2<*<) + 

+*(^) 1/4 (C*X a + X^K^nbexpi-mL-W^d^d^dx^dx 
To calculate the integral (|12]) we make the following change of variables: 

<f% -> exp{-i J D Q n 2 6 2 L- 1 }0^, -> exp{iD a n 2 b 2 L- 1 }^* (13) 
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Then the integral (|P2|) acquires the form 



N 

I = / exp{6 £ E b- 2 K a + i€ + - 2C°C] + (14) 



n=—N a 

2m 



exp { W a b 2 n 2 L -ly x a + x »* exp {-iD a b 2 n 2 L^^) X 

71 Li 

2nbexp{-mL- 1 b 2 n 2 }}d(f)^d(j)yx a *dx a 

Now the quadratic form in the exponent does not depend on D a and therefore the 
corresponding determinant is a trivial constant. So we can calculate the integral 
by finding the stationary point of the exponent, which is defined by the following 
equations: 



i(^) 1/4 X* a 2nbexp{-(tD a + mWn 2 ^ 1 } = 0, n^-N 

&- 2 (0n+l + 2<)+ (15) 

+i(^-) 1/4 X a 2nbexp{(iD a - mWn 2 ^ 1 } = 0, n + -N 
ttL j 

b~ 2 {<j)T N+1 - 2<j>*? N ) + z(^) 1/4 x* Q 2A^6exp{-(2 J D a + m^N 2 ^ 1 } = 



"7TL'- 

2m 
'Id}- 



b-\r-N+i ~ W-n) + i{^) X '\ a 2Nbex V {{iD a - m^N'L- 1 } = 



r- N = <t>- N = o, r N+ i = r N *+i = o 

These equations are most easily solved for small b, when they can be approxi- 
mated by the differential equations: 

<f a + i(^) 1/4 x* a 2texp{-(tD a + m)t 2 L- 1 } = (16) 
'<j) a + i{^j) l/ \ a 2tex V {(iD a - m^L- 1 } = 

r(§) = <n-§) = o, a *(§) = <r(-^) = o 

The solution of these eq.s is 

« = ,t ^S^ #( ' ri, ' i "' +iy)i/ ' ) + < i7) 

( 2 '")' /4 cy. B r L{m + iD ° ] }2tL- 1 - iv» (2 " mL)1/4 



(7TL) 1 /4( m + 2j D«)2 4 2(772 + ^)3/2 

(2mvrL) 1 / 4 
2(m-iD Q ) 3 / 2 ' 



4 



(2m) 1 / 4 -L(m-iD a ) x a (2TxmL) 1 ^ 

+tX (7rL)i/4( m _, jD «)2 ex Pi 4 l 2tL 2{m-iD")W 

where is a Fresnel integral. Substituting this solution into the integral (p~5|) we 

get 

f ^ / >+L / 2 

Integrating by parts and using the fact that </>" t (L/2) = L/2) = 0, we can 
transform this integral to the form 

/ = / exp{£ j L ' 2 dtr s ;r st }d X a *d X a = (20) 

I exp{-E(^) 1/2 f+m^ r, dteM-2mL-H*}+ 
J a 7cL m z + \D a y J-L/2 

A . 2rxnf L ( m+iDa ) } r L/2 

+ 77 2 * , 2 , e^ } + Pi -4- 1 / ^exp{-(m-z J D")L- 1 t 2 }+ 
L(m 2 + (L> a ) 2 ) (m + iD a )L J -l/2 1 v 7 J 

2exp{- L(m 7 Da) } 
(m - iD a )L 
Integrating over \ one gets 



/ dtexp{-(m + iD a )L-H 2 }]} 

J-L/2 



I = Y\(m 2 + (D a ) 2 )(l + 0((mL)- 3 / 2 exp{-^})) (21) 

So if m ~ a -1 , and L = 2iV&, the corrections are of the order 

OiiNba' 1 )- 3 / 2 expi-Nba- 1 }) (22) 



At the same time it is easy to show that replacing the sum in the eq. ( |1^) by the 
integral over t produces corrections of order 0(b 2 a~ 2 ). Therefore taking N ~ 3a6 _1 
we can make the corrections due to the finite size of a lattice less than 1 %. 

Taking into account that 

( l5 D) 2 = -(D) 2 (23) 

we see that 

/ = det(D + m) 2 + 0(b 2 a' 2 ) + 0((mL)- 3 / 2 exp{-^}) (24) 



To conclude we note that our equation @ is a linearized version of the integral (|i~Tp. 
The exponent in the quadratic form of the action (|ll]) may be written in the form 

exp{-i-/ 5 Db 2 L- 1 (2n + 1)}, 2nb < L (25) 

Having in mind that 

II75-DH < 8a" 1 (26) 



one sees that 



l5 Db 2 (2n + l)\\ 



< 8ba- L (27) 

and if 86 < a one can replace the exponential in the action (|H]) by the first terms 
of it's Taylor series, leading to the eq.(|9|). 
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3 Discussion 



We proved that the purely bosonic integral @ is equal to the square of the deter- 
minant of covariant Dirac operator up to the terms 

0(b 2 a- 2 ) + 0((mL)-lexp{-^}) (28) 

Convergence of this approximation is much faster than in the ref. [flj, where the finite 
size effects were supressed polynomially. For finite lattice spacing b the corrections 
decrease exponentially when the number of extra fields N increases. To get the same 
accuracy one needs much smaller number of bosonic fields, which hopefully will 
simplify Monte-Carlo simulations. It is worthwhile to note that to get exponential 
damping it is crucial to introduce the exponential factor into the constraint (|10D. 
As for the choice of the quadratic form in the effective bosonic action, there exists a 
certain freedom. We choose the particular form (O)to simplify analitic calculations. 
Presumably the choice of the quadratic form done in our previous paper [|IJ is also 
possible. We also note that contrary to our previous algorithm [TJ the present 
construction cannot be reduced to Liisher's algorithm with a particular choice of a 
polynomial due to explicit dependence of the constraint on the mass. 
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